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Many-body localization and mobility edge in a disordered Heisenberg spin ladder 
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We examine the interplay of interaction and disorder for a Heisenberg spin ladder system with random fields. 
We identify many-body localized states based on the entanglement entropy scaling, where delocalized and 
localized states have volume and area laws, respectively. We first establish the quantum phase transition at a 
critical random field strength he ~ 8.5 ± 0.5, where all energy eigenstates are localized beyond that value. 
Interestingly, the entanglement entropy and fluctuation of the bipartite magnetization show distinct probability 
distributions which characterize different quantum phases. Furthermore, we show that for weaker h, energy 
eigenstates with higher energy density are delocalized while states at lower energy density are localized. This 
defines a mobility edge and a mobility gap separating these two phases. By following the evolution of low 
energy eigenstates, we observe that the mobility gap grows with increasing the random field strength, which 
drives the system to the phase of the full many-body localization with increasing disorder strength. 

PACS numbers: 75.10.Pq,71.30.+h, 73.22.Gk 


I. INTRODUCTION 


Anderson localization theory' predicts that noninteracting 
electrons are generally localized in one and two-dimensional 
(ID and 2D) disordered systems without either a magnetic 
field or spin-orbit coupling due to destructive quantum in¬ 
terference. It is generally believed that low energy states re¬ 
main localized for weakly interacting systems with char¬ 
acteristic features different from noninteracting systems. Re¬ 
cently, there is renewed interest to examine the Anderson lo¬ 
calization for interacting systems, where the phenomenon of 
many-body localization (MBL)*’® has attracted intense stud¬ 
ies. Many remarkable properties of an MBL phase has been 
established*^® based on combined theoretical and numerical 
studies. For disordered interacting systems, the random disor¬ 
der can drive a dynamic quantum phase transition*’^^’"'^ from 
a delocalized state to an MBL phase, where energy eigen¬ 
states at finite energy density become localized. From the 
quantum information perspective, energy eigenstates in an 
MBL phase have suppressed entanglement entropy satisfy¬ 
ing an area law*’'^’^®’*® scaling with the subsystems boundary 
area in contrast to the volume law scaling expected for an er- 
godic delocalized state. As a consequence, the MBL phase is 
non-ergodic and can not thermalize"’^*’''®, which also chal¬ 
lenges the fundamental “eigenstate thermalization hypothe¬ 
sis” (ETH) for quantum statistical physics®". The MBL state 
may exhibit quantum order or topological order'®’*®’"'®’® at 
finite temperature as excitations at finite energy density are lo¬ 
calized. A phenomenological study'"' further establishes that 
the MBL phase behaves like integrable systems, respecting 
extensive numbers of local conservation laws '5,20,56 
quantum phase transition from an MBL phase to a delocal¬ 
ized ergodic phase may be continuous characterized by a jump 
of the entanglement entropy in the thermodynamic limit'*, 
where both entropy and its variance grow with the system vol¬ 
ume at the critical point'"'’*®. Interestingly, it is conjectured 
that an MBL phase can also have a continuous localization- 
delocalization transition to a new state, where the delocalized 
phase is non-ergodic whose volume law entanglement entropy 
tends to zero as the transition is approached'*. It may be possi¬ 


ble to have the MBL phase in multi-component systems with¬ 
out random disorder®*’®*. The MBL phase may be detected 
experimentally in cold atom systems '2.13,18,19 

So far, much of the quantitative understanding of MBL 
systems are based on numerical exact diagonalization (ED) 
studies'*’'"’*'^*®’®" of spin and electron systems, where the 
dynamic quantum phase transition between a delocalization 
phase to an MBL phase has been demonstrated for differ¬ 
ent ID model systems with spin (or particles) numbers in 
the range of = 10 — 22*®’"'®. There are also some recent 
developments®"^® using tensor network and density matrix 
renormalization group approaches to study such systems. One 
of the conceptually important and unsettled issue is if the mo¬ 
bility edge exists for microscopic system to separate the low 
energy localized state from the higher energy extended states. 
On the one hand, these ED studies*®’"'® have demonstrated the 
energy density dependence of the critical random field, con¬ 
sistent with the existence of the mobility edge. In particular, 
Luitz et al."'® studied the ID Heisenberg chain in a random 
field using the shift-inverted spectral transformation method 
dealing with up to 22 spins, where the finite-size scaling has 
been demonstrated with convincing accuracy supporting the 
existence of the mobility edge. However, the recent numeri¬ 
cal linked cluster expansion study®® for a thermodynamic sys¬ 
tem finds that a higher disorder strength (as the lower bound) 
is required to enter the MBL phase than that obtained by ED 
studies. The reason for such a discrepancy remains not under¬ 
stood. On the theoretical side, it is not clear**’"'"'’®*^" if some 
spatial region with higher energy density may play an impor¬ 
tant role with more extensive wavefunctions, which may melt 
the lower energy eigenstates in the system into delocalized 
states with increasing the system size. Some insight on this is¬ 
sue may come from the earlier study of interacting many-body 
systems with random disorder*"’*', which exhibit the fraction¬ 
alized quantum Hall effect. In such a system, we have demon¬ 
strated that low energy states below a mobility edge have topo¬ 
logical order protected by a mobility gap which separates the 
low energy localized insulating states from the metallic states 
above the mobility gap. These earlier studies suggest that it is 
possible to follow the evolution of the low energy eigenstates 
in disordered interacting systems to detect if the mobility gap 
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generally exists for MBL systems. 

In this paper, we numerically examine the interplay of in¬ 
teraction and random disorder field for two-leg Heisenberg 
ladder systems, which stands between ID and 2D systems^®, 
with the latter being much harder to be systematically stud¬ 
ied based on ED method. We identify MBL states based on 
the bipartite entanglement entropy scaling, and the spectral 
statistics of many-body energy levels. We first establish the 
quantum phase transition at a critical random field strength 
he ^ 8.5 ± 0.5, where all energy eigenstates are localized 
beyond that value. Interestingly, the entanglement entropy 
shows distinct probability distribution in different quantum 
phases, while transition is associated with the divergent vari¬ 
ance for the entropy at the thermodynamic limit^®’^®. Despite 
the small system sizes we can access with the number of spins 
iV = 12 — 20, our entropy distribution matches to the theo¬ 
retic prediction®^ in both delocalized Griffiths phase and MBL 
phase indicating we were able to access universal characteris¬ 
tics of these different quantum phases. Lurthermore, we show 
that at weaker h, energy eigenstates with higher energy den¬ 
sity are delocalized while states at lower density are localized, 
which defines a mobility edge separating these two dynam¬ 
ically distinct quantum phases in agreement with earlier re¬ 
sults for ID spin chain systems^®’'*®. The associated mobility 
gap (the excitation gap of the lowest energy extended state to 
groundstate) may be relevant for experimental measurement 
as it plays the role of the activation gap^®’^^ for quantum trans¬ 
port. By following the evolution of low energy eigenstates, we 
observe that the mobility edge moves to higher energy density 
with the increase of the random field strength, which eventu¬ 
ally drives the system to the phase with full MBL where all 
energy eigenstates are localized. The remaining of the pa¬ 
per is organized as following: In Sec. II, we first introduce 
the two-leg ladder spin model with random fields and briefly 
discuss the method we use to study the system. We present 
the evidence of the MBL phase determined by the entangle¬ 
ment entropy and the energy level statistics studies. In Sec. 
Ill, we study the characteristic features of the quantum phases 
and the transition between the delocalized phase and the MBL 
phase. We also present the evidence for the mobility edge sep¬ 
arating the low energy localized phase from the higher energy 
extended states. Linally, in Sec. IV, we summarize our main 
results and discuss open questions. 


II. SPIN MODEL AND TRANSITION TO AN MBL PHASE 

We study the Heisenberg two leg ladder spin system on the 
square lattice with the following Hamiltonian: 

ihj) 

where the summation {i,j) runs over all distinct nearest 
neighbor bonds with antiferromagnetic coupling J, which is 
set as the units of the energy J = 1. The hi is a ran¬ 
dom magnetic field coupling, which distributes uniformly be¬ 
tween window (— h) with h as the strength of random fields. 


The number of sites of the ladder system can be written as 
N = NxNy with Ny = 2 and Nx is the number of sites along 
each spin chain. 

We perform Lanezos ED calculations to obtain energy 
eigenstates around a fixed value E determined by the tar¬ 
get energy density e for systems with the number of sites 
N = 12 — 20 in the total Sz = 0 sector. Specifically, for 
each disorder configuration, we first calculate the ground state 
energy Eq and the maximum energy which are used to 
define the target energy density e = {E — Eq)/{E max — Eo)- 
We perform more than 1000 disorder configuration average 
for most systems we studied. Physical quantities^® including 
the bipartite entanglement entropy, energy level statistics and 
bipartite fluctuation of the subsystem magnetization are ob¬ 
tained and averaged over disorder configurations and some¬ 
times also averaged over 30 energy eigenstates with energies 




FIG. 1: (Color online) (a) The ratio of entanglement entropy over the 
number of system sites S/N for different systems from V = 6 x 2 to 
9 X 2 at the energy density e = 0.5 as a function of random disorder 
strength h. Curves for different N approximately cross each other 
around a critical random field strength he = 8.5 ± 0.5. (b) The adja¬ 
cent gap ratio r for states with energy density e = 0.5 as a function of 
disorder strength h. Here we see that on small h side, r approaches 
the Gaussian orthogonal ensemble value (0.5307) representing de¬ 
localized states, while at larger h side, r reaches the Poisson value 
(2ln2 — l ~ 0.3863) for larger systems representing localized states. 
All curves cross around the critical value he = 8.0 ± 0.7. 
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closest to the given energy density e as detailed below. 

The bipartite entanglement entropy has been extensively 
used as an effective tool to characterize quantum phases 
for such an interacting system®’^®’"^® We compute the Von 
Neumann entanglement entropy of the ladder system from 
all eigenvalues of the reduced density matrix pA S = 
—Trp^ In Pa, by partitioning the system in the middle along 
the vertical direction (the lengths for two subsystems A and B 
are the integer-parts of Nx/2 and {Nx + l)/2, respectively). 
For an interacting system with weak disorder, the entangle¬ 
ment entropies of higher energy eigenstates are expected to 
follow the volume law and these states are ergodic satisfying 
the ETH^’®. This is in contrast to the behavior of the ground 
state, where the entanglement entropy follows the area law 
(with possibly up to the logarithmic correction depending on 
if there are gapless excitations)'^. By varying the disorder 
strength h, one can detect the possible quantum phase transi- 
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FIG. 2: (Color online) (a) The probability distributions of the bi¬ 
partition entanglement entropy P{S) for spin system at energy den¬ 
sity e = 0.5 with different disorder strengths h = 4, 6, 8 and 10 
for system sizes N — 6 x 2 and 8x2. These results illustrate 
that for stronger disorder case, the distribution always has a long 
tail into higher S values, while for smaller h the long tail is on 
small S side, (b) The probability distributions P{F) of the vari¬ 
ance F =< (SJ)^ > — < Sa (in units of the square of 
the Planck constant h^) of the magnetization of the half system A 
for spin system at energy density e = 0.5 with different disorder 
strengths /i = 4, 6, 8 and 10 for V = 6 x 2 and (b) 8x2. 


tion from the behavior of the entanglement entropy^®’''®. As 
shown in Fig. 1(a), we plot the ratio of entanglement entropy 
over the number of system sites S/N for different systems 
from Nx = 6 to 9 at the energy density e = 0.5 as a func¬ 
tion of random field strength h. On the smaller h side, we see 
the ratio S/N increases with system sizes N and approaches 
a constant indicating the volume law growth of S. With vary¬ 
ing h, all data points approximately cross each other around 
a critical value he ^ 8.5 ± 0.5. On the larger h side, S/N 
approaches zero indicating the low entanglement and non- 
ergodic behavior where energy eigenstates are localized. The 
ladder systems we study here have stronger finite size effect 
(from the even-odd effect of Nx) than the ID spin chain sys¬ 
tems, which is the reason that not all curves cross at the same 
point in Fig. 1(a). 

We further use the level statistics analysis from the ran¬ 
dom matrix theory^'’^^ to probe the localization-delocalization 
characteristics of energy eigenstates. In the delocalized 
regime, the level-spacing distribution is described by the 
Gaussian orthogonal ensemble (GOE) statistics, which rep¬ 
resents extended levels with level-repulsion between them be¬ 
cause of the overlap of energy eigenstates in real space. In the 
localized regime, the level-spacing distribution is determined 
by Poisson statistics as wave-functions close in energy are ex¬ 
ponentially localized with no level repulsion between them^"^. 
In the energy spectrum analysis''®, we define the energy gap 
Sn = En — En-i as the energy difference between the n-th 
and (n — l)-th eigenstates, then the adjacent gap ratio can be 
defined as r„ = min(Sn,Sn+i)/'rnax(dn,Sn+i)- We average 
the gap ratio r =< > over states near the spectrum cen¬ 

ter at e = 0.5 for 30 eigenstates and 1000 random disorder 
configurations for each given disorder strength h. As shown 
in Fig. 1(b), we see that at the small h side, r approaches 
the Gaussian orthogonal ensemble value (0.5307) represent¬ 
ing delocalized states, while at stronger h side, r reaches the 
Poisson value (2/n2 — 1 ~ 0.3863) for larger systems rep- 
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FIG. 3: (Color online) (a) The variance (AS)^ =< > — < 

S of the entanglement entropy at energy density e = 0.5 for dif¬ 
ferent h. AS reaches the peak value at hp smaller than the identified 
he for quantum phase transition. Clearly, hp shifts to higher h with 
the increase of N. 
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resenting the level statistics of localized states. All curves 
cross around the critical value he ^ 8.0 — 8.5. Due to the 
stronger finite size effect coming from the even-odd effect for 
Nx, we are not attempting to do a finite-size scaling. Instead, 
we will focus on the universal behavior of the different phases 
to explore the nature of the quantum phases and the transition 
involved here. 
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FIG. 4; (Color online) The entropy of each energy eigenstate Si 
for low energy eigenstates averaged over disorder configurations is 
shown as a function of its average energy density. Si/N curves for 
different system sizes A = 6x2, 8x2 and 10 x 2 cross at a criti¬ 
cal energy density Sc, which separates the higher energy delocalized 
states from the lower energy localized states. 


III. NATURE OF QUANTUM PHASES AND PHASE 
TRANSITION 

A. Probability distributions of entanglement entropy and 
variance of bipartite magnetization 

Now we turn to the study of the probability distribution of 
the entropy P{S) for spin system at energy density e = 0.5 
for different disorder strengths ft, = 4, 6, 8 and 10 crossing 
two different phases obtained for ensembles with 30 energy 
eigenstates and 1000 disorder configurations for each ft. As 
shown in the Fig. 2(a), on the small ft side, we see that the 
peak position of the distribution P{S) (which reflects the av¬ 
erage of S) moves to the larger S value with increasing sys¬ 
tem size N, indicating a consistency with the volume law for 
the entropy. Close to the transition point for ft = 6 we find 
that the distribution P{S) becomes much broadened while the 
peak position moves to the higher S with the increase of N, 
but the peak height reduces at the same time. As we move to 
the higher ft side, we see that the distribution again becomes 
sharper, with two peaks showing for each P{S) curve which 
may be related to the non-ergodic character of the localized 
phase. Furthermore, we also see that for the stronger disor¬ 
der case, the distribution always has a long tail into higher S 
values, while for smaller ft the long tail is at the smaller S 
side. To compare with recent theoretical description of the 
MBL of ID system®^, we find that the entropy distribution 
is very similar to the ones obtained based on their real space 
renormalization group simulations. Specifically, at ft = 4 our 
distribution has a long tail in small S region, which matches 
to the one for the Griffiths phase just before the transition to 
the MBL phase. At ft = 10, we find the P{S) is peaked at 
S' = 0 and shows an exponential decay tail on the larger S 
side. 

We compare the entanglement entropy behavior with the bi¬ 
partite fluctuations F of the subsystem magnetization 
which is defined as F = ((S;|)^) — as shown in Fig. 

2(b). The distribution P{F) exhibits very similar behavior as 
P{S) for ft closer to the quantum phase transition and in the 
MBL phase. The similar two-peaks structure is also clear for 
P{F) on the larger ft side in the MBL phase. The only differ¬ 
ence worth mentioning is that with weaker disorder ft = 4, the 
P{F) demonstrates the normal Gaussian distribution, which 
is sharp and near symmetric about the peak. 

The variance of the entanglement entropy has been shown 
to be an excellent quantity^®’®^ for identifying the quantum 
phase transition from ID spin chain studies. Here we show the 
variance AS' of the entanglement entropy averaged over 30 
different energy eigenstates around the energy density £ = 0.5 
and 1000 disorder configurations. In agreement with these 
observations for ID systems^®’"^®, we find that the AS is small 
in both small ft and large ft sides, and demonstrates a peak for 
the intermediate ft as shown in Fig. 3. We observe that the 
peak value of AS increases with N, which may diverge at the 
transition point. The position of the peak hp is smaller than the 
previously identified and it shows a trend of approaching 
he with the increase of N. These results are consistent with 
the phenomenological theory®^ established based on the real 
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space renormalization group studies, which indicate that we 
are observing intrinsic properties of the MBL phase and the 
related quantum phase transition for system sizes we study. 

B. Mobility edge and mobility gap 

We have shown that the disorder can drive a quantum phase 
transition, where all states near the center of energy spectrum 
are localized. In fact, all other states with different energy 
density are also localized, thus we enter the full MBL phase 
tuned by h (see Fig. 4(c) as an example). To address the 
issue if the mobility edge naturally exists in such a system 
separating low energy states with the area-law entanglement 
entropy from higher energy states with volume-law behavior, 
we follow lower energy eigenstate by obtaining hundreds of 
these eigenstates in ED. We follow the entropy of each energy 
eigenstate Si in the lower energy density regime and average 
that over 1000 disorder configurations. 

As shown in Fig. 4(a-b) for h = 2 and 4, we iden¬ 
tify that the disorder configuration averaged entropy Si for 
the i — th energy eigenstate is a smooth increasing func¬ 
tion of eigen energy Ei or its average energy density Ei =< 
{Ei - Eq)/{E^ ax - Eq) >. The entropy per site Si/N for 
different system sizes iV = 6x2, 8x2 and 10 x 2 crosses 
around a critical energy density Ec, which separates higher en¬ 
ergy states with volume-law entropy (delocalized states) from 
lower energy localized states with Si/N approaching zero vi¬ 
olating the volume law. Here all the data we show have even 
Nx = 6, 8 and 10 with reduced finite size effect. The cross¬ 
ing point determines the mobility edge. With increasing h, 
the entropy of the low-lying eigenstate is reduced and the mo¬ 
bility edge is being pushed to the higher energy density from 


Ec 0.075 aX h = 2 to Ec 0.1 ath = A. As shown in Fig. 
4(c), we further move to the stronger disorder case, ath = 10, 
we see that S/N at different energy density (here we averaged 
over both disorder configurations and energy eigenstates for 
each energy density e) is always a decreasing function with 
increasing N demonstrating all states are localized. 


IV. SUMMARY AND DISCUSSION 

We have identified the disorder driven dynamic quantum 
phase transition from ergodic delocalized phase to an MBL 
non-ergodic phase for two-leg ladder Heisenberg spin systems 
with random field disorder. The characteristic distributions 
of the entanglement entropy for both the delocalized Griffiths 
phase and the MBL phase agree with the theoretical descrip¬ 
tion for the MBL®^. Furthermore, we show that for weaker h, 
energy eigenstates with higher energy density are delocalized, 
while states at lower density are localized, which defines a 
mobility edge separating these two dynamically distinct quan¬ 
tum states in agreement with earlier results for ID spin chain 
systems^®’'^®. On the quantitative side, we find that the Heisen¬ 
berg ladder requires a much higher critical disorder strength 
compared to the ID spin chain model^® to enter the full MBL 
phase. It is highly desired to further study multi-leg ladders 
and explore scaling behaviors to the 2D limit, which we leave 
for a future study. 
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